
Calhoun 

Iniiiiuiiortfl Arthivcof (he Navjl Pwigndualt School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 



Theses and Dissertations 


1. Thesis and Dissertation Collection, all items 


2012-12 

Design of a maximum power point tracker 

with simulation, analysis, and comparison of algorithms 


Wurz, Jeff 

Monterey, California. Naval Postgraduate School 


http://hdl.handle.net/10945/27924 


Downloaded from NPS Archive: Calhoun 



DUDLEY 

KNOX 

LIBRARY 


htt p ://w w w. n ps.e-du/l ib ra ry 


Calhoun is the Naval Postgraduate School's public access digitaI repository for 
research materials and institutional publications created by the NPS community. 
Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 
appointed —and published —scholarly author. 

Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 



NAVAL 

POSTGRADUATE 

SCHOOL 


MONTEREY, CALIFORNIA 

THESIS 


DESIGN OF A MAXIMUM POWER POINT TRACKER 
WITH SIMULATION, ANALYSIS, AND COMPARISON OF 


ALGORITHMS 


by 


Jeff Wurz 


December 2012 

Thesis Advisor: 

Robert Ashton 

Second Reader: 

Vladimir Dobrokhodov 


Approved for public release; distribution is unlimited 




THIS PAGE INTENTIONALLY LEFT BLANK 









THIS PAGE INTENTIONALLY LEFT BLANK 


11 



Approved for public release; distribution is unlimited 


DESIGN OF A MAXIMUM POWER POINT TRACKER WITH SIMULATION, 
ANALYSIS, AND COMPARISON OF ALGORITHMS 


Jeff Wurz 

Civilian, United States Navy 

B.S., Electrical Engineering, University of California, Davis, 2008 

Submitted in partial fulfillment of the 
requirements for the degree of 

MASTER OF SCIENCE IN ELECTRICAL ENGINEERING 

from the 

NAVAL POSTGRADUATE SCHOOL 
December 2012 


Author: 


Jeff Wurz 


Approved by: Robert Ashton 

Thesis Advisor 


Vladimir Dobrokhodov 
Second Reader 


Clark Robertson 

Chair, Department of Electrical Engineering 



THIS PAGE INTENTIONALLY LEFT BLANK 


IV 



ABSTRACT 


In this thesis, the advantages of three different maximum power point tracking ( MPPT ) 
algorithm are investigated. By simulation, the performance and efficiency of these algo¬ 
rithms was analyzed. By using MATLAB’s SimPowerSystems block set, we created the 
model comprised of a Kyocera KD135GX-LP solar panel powering a buck converter con¬ 
trolled by the MPPT algorithms driving a resistive load. The main objective was to track 
the maximum power point ( MPP ) of the solar array by modulating the buck converter’s 
duty cycle, thereby, optimizing the power output of the panel. The three algorithms ob¬ 
served performance was on par with other real world tests of these algorithms as seen in 
other published work. The Perturb and Observe ( P&O ) algorithm performed with a higher 
overall efficiency and was able to track the MPP quickly, while the Incremental Conduc¬ 
tance (InC) algorithm had similar performance but requires more intensive calculations. 
The analysis of these algorithms led to a greater understanding of where the inefficiencies 
of this type of system are located, allowing improvement in future work on this subject. 
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Executive Summary 


The objective of this research was to investigate the improvements in solar power genera¬ 
tion through use of a maximum power point tracker ( MPPT ) in conjunction with a buck 
converter. Solar energy is all around us during the day, the only problem is the inefficiency 
of current solar cells. Every year advances are made in the efficiency of solar cells, but 
overall the cells’ greatest source of correctable inefficiency is in the implementation of the 
MPPT. 
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Figure 1: This is the Current versus Voltage and Power versus Voltage plots from the 
Kyocera KD135GX-LP panel modeled. 


Solar cells operate as a reverse biased diode. Solar cells absorb light incident on their 
surface which recombines in the active region of the diode. In order to extract the most 
power from a solar cell, the equivalent resistance which allows the most voltage and current 






















output from the cell must be maintained. The maximum power point ( MPP ) can vary 
depending on the magnitude of irradiance and the temperature of the solar cell as seen in 
Figure 1. The major differences in the algorithms which will search for the MPP are speed, 
overshoot, and tracking efficiency. These three parameters must be balanced in order to 
achieve optimal power transfer within the system. 

The optimized solar cell output power flows into a bucking converter which has its set point 
determined by the MPPT algorithm. There are three specific MPPT algorithms which were 
investigated: Perturb and Observe ( P&O ), Incremental Conductance (InC) and Constant 
Duty Cycle. By analyzing the output of these algorithms under different known conditions, 
the performance of each of these controllers was assessed. After generating the code nec¬ 
essary to implement these algorithms in the MATLAB environment, it was necessary to 
also construct an appropriate model using the SimPowerSystems block set to simulate the 
system, as seen in Figure 2. 



Figure 2: This figure shows the functional portion of the model under test in Simulink, 
which includes the solar array, buck converter, and resistive load. 

The major differences in the algorithms which search for the MPP are speed, overshoot, 
and tracking efficiency. These three parameters must be balanced in order to achieve an 
optimal system. 
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The control of a bucking converter is achieved by modulating the duration the switch is 
‘on’ or ‘off’ as seen in Figure 3. Essentially, the MOSFET switch controls the power flow, 
while the output low pass LC filter removes the switch discontinuity. Through control of 
the duty cycle (D). the output voltage is equal to the input voltage divided by I). Tuning 
the output filter reduces the output voltage ripple to an acceptable level. 



Figure 3: Bucking DC-DC converter showing simple circuit diagram (a), output voltage 
and frequency spectrum of output (b). 
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The comparison of the MPPT algorithms was a worthwhile endeavor. This simulation 
allowed for optimizations to be made prior to purchasing components or building a physical 
circuit. Pending future work, this effort may have great impact on solar powered unmanned 
air vehicles ( UAV ). 

By optimizing the buck converter independent of the MPPT and the solar array, the effi¬ 
ciency of the converter was greater than 94% during steady state operation. This makes the 
conversion from the solar array quite efficient and attractive. 

The ability to track the MPP further optimized the low efficiency solar array output. Since 
the solar array is the source of power to the system, any improvement in the ability to track 
the MPP similarly improves the overall system performance. 

The actual data gathering runs were done in three parts. The first part included the P&O 
algorithm with constant 1000 W j nr Irradiance, a step in irradiance from 600 W/nr to 
1000 W/m 1 at 0.08 5 to 0.01 s, and a random changing of irradiance to the solar array. 
Next, the constant D controller was used as well with the same three input signals. The 
constant D was used as the reference for comparison with the other algorithms. Finally, 
the InC algorithm was tested; although, the tuning of the parameters for this algorithm 
was the most difficult. The algorithm leads to either a very fast response to transients or 
undershoots the MPP. Because of the difficulty optimizing and the amount of time required 
to actually run each test, the InC algorithm was not optimized to match the performance of 
the P&O algorithm. 

An innovation in the P&O algorithm was to include a two stage change in D in order to 
converge quickly to the target, yet still have fine granularity when approaching the MPP. 
This was not very difficult to implement and overall provided very fast tracking and con¬ 
vergence to the MPP. In the Appendix, the code used for this algorithm is provided. It 
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can be seen that with a multistage MPPT , the overall speed of the system remains fast and 
tracks the MPP well. Since this is not division, the operation does not require many extra 
calculation cycles to complete. 

Overall, the expected outcome of the P&O algorithm was quite good. Since the algorithm 
can be optimized to meet certain timing requirements, it responded faster and performed 
better than the InC algorithm. This can only be accounted for due to the complexities of the 
InC algorithm computations and speed of operation. Due to the difficulty in optimizing the 
InC algorithm, the results from those runs do not accurately tell us the optimal operation 
of the algorithm. Given more time to properly optimize AD, the algorithm should have 
tracked the MPP quicker and with more certainty than the P&O algorithm. 
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CHAPTER 1: 
INTRODUCTION 


It is a warm summer day. You feel the sun warm your skin and rejuvenate your motiva¬ 
tion. The sun generates more energy than anything else in the solar system, but due to our 
distance from this writhing fireball, we are perfectly comfortable and able to sustain life. 
Life has figured out a way to harness this power emanating from 93 million miles away. By 
harnessing this power that is available all around us, we enable ourselves to free our bodies 
from the toils we once endured not 40 years ago. More and more automation is taking over 
our lives, and to power that automation we require a source of energy that will not be going 
away anytime soon, the perfect example of this being our sun. Without the sun, life would 
cease to be and, therefore, why not harness the sun to create a better life for those around 
to live it? 


1.1 Objectives 

The objective of this thesis is to determine the most cost effective and efficient method for 
maximum power point tracking ( MPPT ). Many optimizations can be done to a system. 
Knowing which will bring about the greatest change is truly the challenge. 

Since fossil fuels are detrimental to the environment and also non-renewable, there has 
been an upsurge of interest in clean and renewable energy. While more than one option is 
available to fill that void, the most interesting and widespread so far is photovoltaics. Pho- 
tovoltaics are semiconductor devices which convert solar irradiation in the visible spectrum 
to generate direct current (DC). With recent advances in technology and discoveries of new 
materials, solar cells are increasing in efficiency and flexibility. This change in what solar 
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panels are composed of has greatly increased the applications of these panels. The surface 
they are supported on no longer needs to be rigid and planar, and due to this they are now 
able to be used in applications that previously could not use photovoltaics. 

In order to fully harness the power of the sun, we need to optimize the output from the solar 
panels. In the simplest terms a solar panel is a current source when short circuited and has 
a voltage while in open circuit configuration. The key to attaining the greatest output from 
the panels is to maintain an output impedance so that the solar panels operate at a point that 
corresponds to its maximum power point ( MPP ). The MPP is the operating point required 
to optimize the output power of the panels. Depending on weather conditions, the output 
power is proportional to the amount of incident light on the surface of the photovoltaics, as 
well as being inversely proportional to the temperature of the panels. This means that the 
hotter the panels get the less current they can provide at the output. 

In order to convert the power from the solar array, we require a power converter that will be 
able to extract the most power from the array. By connecting a buck converter to the output 
of the solar array, we are then able to control the voltage of the solar array by varying the 
duty cycle (D) of the buck converter. When one is changing D, we want the voltage and 
current to provide the most power at a specific voltage level. In order to reach this voltage 
level, a control algorithm is implemented to track and follow the highest input power from 
the solar array. When this algorithm is functioning correctly, it is said to be an MPPT. 

1.2 Motivation 

Clean and renewable energy has greatly increased funding for supporting areas in research 
and development. Today, energy is necessary for nearly everything. We need a more 
efficient method of transferral and transformation of DC for different voltages and power 
levels. By analyzing the power transfer in a buck converter, we will be able to better 
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understand how and why these converters function so well. 


In our modern society, it has become increasingly necessary to shift away from non¬ 
renewable resources. A change might occur if energy prices were to increase drastically, 
but by then it might be too late to switch our focus and essentially rebuild our machines to 
run off of fuel that is readily available. The effect can be seen today in the car industry as 
currently it is driven by petroleum. With every increase in petroleum prices, a public out¬ 
cry is heard and more and more people are replacing their inefficient vehicles with modern 
hybrid or electric drive vehicles. We can reduce the effect of petroleum pollution by shift¬ 
ing to renewable resources, such as photovoltaics. In addition to this need for renewable 
resources, there is also a similar need for efficient transformation of power from one form 
to another. 

Since solar panels are quite inefficient to begin with, any method of improving power trans¬ 
fer from these panels is a worthwhile gain in overall efficiency. One method to accomplish 
this is to operate the panels at their MPP. By connecting a solar array to a buck converter, 
we can control the output voltage of the panels. The parameter to be controlled is the 
percentage of time a switch is active during each control cycle. In order to compare the 
performance of different algorithms designed to control the buck converter, a simulation 
was made to compare the algorithms. 

1.3 Technology Overview 

In order for the reader to get up to speed, we need to cover the relevant histories of each 
component. First, the history and important discoveries in photovoltaics is included. These 
milestones in solar cell history are important because this will further explain the trend of 
photovoltaic research. Then we will discuss the buck switch mode power converter and its 
operation. Finally, the chapters’ contents will be described. 
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1.3.1 Solar Cells 


The late 1800s were the beginning of man’s search for the electrical properties of materials. 
Eleven years after selenium was discovered to have photoconductivity in 1873, the first 
solar cell was patented by Edward Weston under patent number US389125 [1]. He was 
only the first in a long line of scientists who were interested in photoconductivity as well 
as the photovoltaic effect. Albert Einstein published a paper explaining through quantum 
physics the photoelectic effect [2]. His paper helped many other scientists to understand 
the mechanism of photovoltaic action and, therefore, this technology has improved over 
time. Bell labs was the first to lead industry by developing a silicon P-N junction photocell. 
These first photocells have an efficiency of six percent [3]. 

The first solar powered satellite was the TIROS-1 launched in 1960 [4]. Previous satellites 
used batteries that were charged before heading to space and had no way to be recharged in 
space. Since technology for batteries was in its early stages of development, the power ca¬ 
pacity for these small communications relays was quite low until they had a way to recharge 
the batteries while in space. Seven years later, in 1967, Soyuz 1 was the first manned space¬ 
craft to have its power replenished by solar panels [5]. From that point onward, nearly all 
satellites have been solar powered. 

In 1980, thin film solar cells were developed by a team at the University of Delaware [6]. 
Thin film cells can use one of four different absorbing semiconductors: amorphous — Si, 
CdTe, CuInSe 2 , and CiijS. Depending on which absorbing material is chosen, the effi¬ 
ciency and cost can be greatly affected. The lower cost of manufacturing and low mass 
encouraged more research into thin film solar cells. In many of today’s integrated applica¬ 
tions, thin film solar cells are an important component. Their importance is derived from 
their flexibility and their thinness. Current bulk-silicon-based photovoltaics have a very 
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thick profile and cannot be bent or deformed even a small amount without shattering. 


Recently, there have been many innovations that have allowed photovoltaics with efficien¬ 
cies greater than 40% to be created by the University of Delaware [7] . While these are still 
research cells, it says a lot about the future of solar since they have already improved these 
efficiencies from 10% to 40% in 50 years, as can be seen in Figure 1.1. The time needed to 
start a new fabrication line as well as the research needed to further optimize new types of 
manufacturing is usually around ten years; this assumes that production costs are reason¬ 
able and their special structures used to achieve the increased efficiency are reproducible. 


Best Research-Cell Efficiencies 
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Figure 1.1: Best research solar cells in production by year. Image from [8]. 


When analyzing a solar cell, efficiency is one of the greatest parameters which is measured. 
These measurements are done with special lights, solar simulators, which will illuminate 
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the array with a constant 1000 W/nr as we see on the surface of Earth. By knowing the 
total amount of power in the light that is incident on the surface of the cells, the efficiency 
can be found from 


ri = — x 100 %, ( 1 . 1 ) 

Pin 

where P out is the output power and P m is the input power. 

In order to get the correct efficiency measurement of the solar array, we need to normalize 
based on area. By dividing energy by the area, we obtain the normalized energy output 


where £ is solar irradiance energy and A is area of the solar array. 

1.3.2 Buck Converter 

Output Voltage Ripple 

It should be noted that in switch-mode DC power supplies, one of the more common metrics 
to judge the performance of the converter is output voltage ripple. Usually, voltage ripple 
is on the order of <1% of the total output voltage. Therefore, analysis assuming v a (t) = V a 
is valid. It should be shown that the output ripple is consistent with the discussion of the 
low-pass filter characteristics of the input [9]. 


1.4 Thesis Organization 

In Chapter 1, background information about solar technology was provided. The current 
problem with power output from solar panels was described. This information was used to 
further solve our primary research questions. An overview of the solar cell was discussed 
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as well as the necessity to convert power between different voltages and currents using a 
buck switch mode power converter. 

In Chapter 2, the literature review of other papers written relating to the following topics 
is presented: buck converters, solar arrays, or MPPT. The differences between this project 
and referenced papers are discussed. Any improvements made to the circuit or simulation 
due to the literature review are noted. 

In Chapter 3, the functionality of each component in the system is explained. They are 
modeled in MATLAB to simulate the real-world components. The operation of each circuit 
component in the system is covered in Section 3.1. The modeling of these components in 
MATLAB is covered in Section 3.2. 

In Chapter 4, test runs on the simulation of the buck converter with MPPT control are used 
to show the efficiency of power conversion. The I-V and P-V curves of the solar panel are 
also illustrated. How the model is assembled in MATLAB along with any challenges along 
the way are also described. It was also necessary to test the performance of the model. 

Chapter 5 contains a comparison of the three major MPPT algorithms and a discussion of 
how each has its own optimal parameters which can be optimized. The effect of tuning the 
MPPT algorithms through both sample rate and gain of duty cycle change are also shown. 
We will also compare the output response from each MPPT algorithm under each of the 
three illumination profiles. 

Chapter 6 is the conclusion of the thesis. Useful ideas to examine in the future are covered 
as well as other uses for this type of power converter. Modifications that can be done in 
order to compensate for changes in design parameters, as well as recommendations to ease 
future work are also discussed. 
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CHAPTER 2: 

LITERATURE REVIEW 


2.1 Masters Thesis - Hurd 

In this thesis, William R. Hurd documented his best effort to use solar cells to increase 
the endurance of a unmanned air vehicle ( UAV ). In order to optimize the output of his 
solar panels, he needed to use an MPPT. The MPPT he chose also contained a battery 
charging circuit. This simplified his designs to allow him to focus on the actual mounting 
and positioning of the solar cells. The panels he used were removed from the manufac¬ 
turer’s encapsulation. He then laminated the cells to the top of the air vehicle’s wing. A 
major difference between his project and the one described in this thesis is the voltages of 
his solar array input and the voltage of his battery system. While his system had a higher 
battery voltage, the one considered in this thesis has a higher solar array voltage. His de¬ 
sign required a boost converter, while the one considered in this thesis requires a buck 
converter. [10] 

2.2 Masters Thesis - Coba 

This thesis was very similar in content to William R. Hurd’s thesis. Javer V. Coba was 
attempting to continue the work of Mr. Hurd. Javier Coba was focused on a specific 
production UAV, the Raven RQ-1 IB. While Mr.Hurd did have the Raven UAV in mind, he 
did not have the chance to test his addition in the field. By their calculations, they nearly 
doubled the endurance of the unmodified UAV. This improvement did come at a cost as the 
solar wing they fabricated was too large to fold and store in the original manned portable 
UAV package. Because of this, they concluded that any solar array extension of the Raven 
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RQ-11B wing system should be done by the company when the wings are manufactured 
so as to provide a reliable and scaled form factor and still enable a two-man crew to deploy 
the UAV in field. [11] 

2.3 Modeling and Simulation of Photovoltaic Module Us¬ 
ing MATLAB/Simulink 

The steps taken to measure the output of two specific solar array modules, MSX60 and 
MSX64, were described in this paper. Measurements were taken in order to impliment a 
model in MATLAB. Datasheet parameters were used as input into the equation describ¬ 
ing the current output of the solar panel. They used the model in order to simulate the 
MPP and graphed the current versus voltage and power versus voltage, showing the higher 
output panel producing more peak power than the lower rated panel. Their experiements 
confirmed that the simulation was well matched to the datasheet information. [12] 

2.4 Advantage of Boost Versus Buck Topology for Maxi¬ 
mum Power Point Tracker in Photovoltaic Systems 

Since the buck converter was used to step down the voltage of the solar array, the drawbacks 
of using buck with a solar array was analyzed. In this paper, they discussed how the use of a 
buck converter is not optimal for solar arrays. The most important aspect of the experiment 
is getting as much power out of the panels as possible. While a buck converter is connected 
to a solar array, it operates the solar array in a state of source discontinuity. This means that 
while the switch is ‘off’, the panel is disconnected and no current is flowing from the panel 
to the converter. This state can be considered pure loss. Without a continuous connection to 
harness the output of the panel, all energy that would have been collected while the switch 
is in the ‘off’ position is lost. The way to circumvent this problem is to add a large input 
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capacitor, which can store at least the same amount of charge that is generated while the 
panel is disconnected from the converter. 

In general the buck converter has a control advantage. While the boost has a right-half 
plane zero, the buck does not. This right-half plane zero causes the response of the boost 
converter to first fluctuate in the opposite direction prior to correcting itself. [13] 

2.5 Comparison of MPPT Algorithms for DC — DC Con¬ 
verters Based PV Systems 

The differences between two popular MPPT algorithms, P&O and InC were described in 
this paper. In order to get a good metric for comparison, the authors designed both power 
conversion circuits to be optimal for that specific algorithm. This was done to compare the 
algorithms on an even footing. The conclusion they draw is that the best controller is InC, 
but if you consider increased cost of implementation of the optimal circuit, then P&O gains 
favor. Due to the lower cost of measurement devices needed for the algorithm. Perturb and 
Observe is the most cost effective and efficient MPPT algorithm. [14] 

2.6 High Efficiency Switched Capacitor Buck-Boost Con¬ 
verter for PV Application 

A switching capacitor buck-boost hybrid converter for module level distributed -MPPT PV 
applications was introduced and described in this paper. The operation principle of the 
proposed converter is covered including the detailed operating of a resonant charge pump 
converter. This converter achieved 92.5% efficiency in the experiment. While not very 
similar to the application described in this thesis, it is good to have a few sources which 
discuss other methods of implimentation of power converters. [15] 
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2.7 A Novel MPPT Charge Regulator for a Photovoltaic 


Stand-alone Telecommunication System 

A commercial MPPT , SIELTE S4007 with their own converter which uses a two-phase 
synchronous buck topology was discussed in this paper. They found that there was a sig¬ 
nificant loss in power when using the SIELTE module to track the MPP. Its error from 
ideal was always greater than 5%. It is interesting to note the way that they implimented 
the two-phase synchronous buck topology. They used Kyocera 125-G2 panel as the basis 
for their circuit simulator model. [16] 

2.8 Improved Circuit Model of Photovoltaic Array 

In this paper, the authors first explain the current circuit model of how a photovoltaic array 
functions. They improve upon the model by using a next order piecewise linear mathe¬ 
matical model. They were able to accomplish this by measuring the output of a real panel 
under known insolation and temperature. By using those measurements as data points on a 
piecewise function to model the photovoltaic array, they showed that this model was able 
to account for mismatched panels connected together as well as partial shading of sections 
of the array. [17] 

2.9 Comparitive Study of Maximum Power Point Track¬ 
ing Algorithms 

The most common algorithm for MPPT is P&O, yet it is explained in this paper that it is not 
the most efficient algorithm. The most common reason for using P&O is the simplicity of 
the circuit. This leads to the negation of need for a complicated processing unit. Their main 
goal was to evenly compare the algorithms across equivalently optimized hardware as well 
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as with the algorithms optimized for that hardware. They conclude that the increase in cost 
of implementation of InC greatly outweighs any improvement in efficiency gained. [18] 


2.10 Design and Implimentation of Maximum Power Point 
Tracking Algorithm for a Standalone PV System 

Using a model based on Simulink blocks in MATLAB, the authors simulated the action 
of an MPPT system. With a variable irradiance from the solar array, the authors use an 
MPPT algorithm to control a boost converter driving a resistive load. In this paper, the 
authors compare P&O with InC. They note that tuning the algorithm using incremental 
changes is important for both algorithms. InC uses far more complicated calculations and, 
thus, requires more powerful hardware when physically implemented. [19] 
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CHAPTER 3: 

SYSTEM OVERVIEW 


3.1 System Diagram and Description 

In this chapter, the use of MATLAB with Simscape, SimPowerSystems and SimElectronics 
toolboxes to simulate the complete system as seen in Figure 3.1 are described. 

Initially, sunlight hits the solar panel array. After being recombined in the active region of 
the solar cell, the electrons flow into a buck converter. The buck converter then lowers the 
voltage and increases the current. Control of this converter is provided by an MPPT algo¬ 
rithm, which monitors the solar array output and uses the current and voltage measurements 
from the array to track the highest power output of the array. 

One such MPPT algorithm will perturb the output setting and then observe by measurement 
the change in output power. When the power reaches the MPP, the algorithm will stop 
changing the setpoint for the duty cycle ( D ). One of the parameters to investigate how it 
effects the overall system is the sampling rate and how often the MPPT will change the 
D of the buck converter. If the sampling rate is too fast, the new D might oscillate around 
the actual MPP , thereby, never extracting the maximum amount of power. On top of that, 
there might be a transient instability created by not allowing the step response of the output 
voltage to settle after changing the setpoint of D. 
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Figure 3.1: The functional portion of the model under test in Simulink, which includes the 
solar array, buck converter, and resistive load. 

The algorithm will use the MPP to set the D of the pulse width modulation ( PWM ) for the 
buck converter. Under proper conditions, the greatest amount of power is transferred to the 
load with minimal losses. The losses of efficiency in the system are due to MPP tracking 
error, buck converter efficiency, and solar cell conversion efficiency. Due to the lack of 
discrete resistive circuit elements, a buck converter can operate with efficiencies greater 
than 95% [9]. The total output power of the buck converter delivered to the load should 
be at least 95% of the input power from the panels. When the panels react to changes in 
irradiance, there is a disturbance in the setpoint needed to keep the panels at their MPP. 
During the time while the algorithm is tracking the new MPP, there is a small amount of 
inefficiency attributed to tracking error. This error is based on the irradiance, making it a 
difficult error to account for in the physical world. In simulation it can be quite easy to 
measure since we can control the illumination of the panel in a controlled manner. The 
tracking error will be tested in simulation to use as a metric for comparison of the MPPT 
algorithms. 

Measurement of the input and output power are the best metric to determine how well 
the buck converter, MOSFET, gate drive, and MPPT algorithm work together. From this 
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measurement of input and output power, efficiency can be calculated using Equation 1.1. 


3.2 Modeling of Components Under Test 

In order to simulate the analog components which are in use for this thesis, we need to 
use mathematical approximations to what the elements actually do. Without this step, it is 
difficult to compare these circuits. By using these models we are then able to do discrete 
analysis with the assistance of a computer for calculations. Calculating these by hand is 
time prohibitive. In MATLAB there are many blocksets already created to ease the use 
of their software environment. By connecting these mathematical models, we are able to 
simulate the circuit and extract close to real-world data. Once the correct tuning has been 
done in simulation, the data can then be use the data to instantiate a real-world device. 

3.2.1 Solar Array 

For this component, we created a model of the solar array based on 

Id = Isat(e&- 1), (3.1) 

where Id is the diode current (A), V)/ is the voltage across the diode (V), I sat is the diode 
reverse saturation current (A), q is the electron charge (1.602 x 10 l9 C), k is Boltzmann’s 
constant (1.381 x 10 ~ 23 J/K), T is the junction temperature (K), and e is Napier’s Constant 
(approximately 2.71828). 

As stated by Mohan: 

The cell characteristic at a given irradiance and temperature basically consists 
of two segments: the constant-voltage segment and the constant-current seg¬ 
ment. The current is limited as the cell is short-circuited. The maximum power 
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condition occurs at the knee of the characteristic where the two segments meet. 
It is desirable to operate at the MPP. Ideally, a pure DC current should be 
drawn from the solar array, though the reduction in delivered power is not very 
large even in the presence of a fair amount of ripple current. To ensure that the 
array continues operating at the MPP , a P&O method is used where at regular 
intervals the amount of current drawn is perturbed and the resulting power out¬ 
put is observed. If an increased current results in a higher power, it is further 
increased until power output begins to decline. On the other hand, if an in¬ 
crease in current results in less power than before, then the current is increased 
until the power output stops increasing and begins to go down. [9] 


For this simulation, a solar array that could output greater than 100 W of power was se¬ 
lected. The panel is made by Kyocera, KD135GX-LP. The datasheet for this panel can 
be found in the Appendix. This panel has a open circuit voltage V oc of 22 V with a short 
circuit current I sc of 8.4A. The voltage at maximum power V mp is 17.7 V while the current 
at maximum power I mp is 7.6 A. By plugging these values into 

P = IxV , (3.2) 

we attain 134.5 W of output power from the solar array under full solar illumination. 

In order to model the varying amount of shading and cloud cover, a profile of illumination 
was devised that will check the MPPT' s ability to change setpoints within a reasonable 
amount of time. 

By varying the incident illumination of the solar array as well as the output resistance, we 
can extract the Current versus Voltage and Power versus Voltage graphs to find the MPP of 
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the solar array, seen in Figure 3.2. The peak of the Power versus Voltage graph is the MPP. 


The mathematical model of a diode, given in Equation 3.1, was used and integrated with a 
controlled current source from SimPowerSystems. With this block, the buck converter can 
modulate the output power. In Figure 3.3, we can see the Simulink diode model used to 
generate the current in the solar array model shown in Figure 3.4. 




Figure 3.2: The Current versus Voltage and Power versus Voltage plots from the Kyocera 
KD135GX-FP panel modeled. 

3.2.2 Buck Converter 

For the buck converter, we modified the standard topology by moving the MOSFET to the 
ground leg of the circuit. This was done to reduce the large voltage spikes which might 
cause damage to the components. The improved circuit is illustrated in Figure 3.5. The 
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Figure 3.3: The model of a diode which is used within the solar array model. 



Figure 3.4: The solar array model used to simulate the Kyocera KD135GX-LP panel. 

input is coming from the solar array. The output is being modeled as a resistive load for 
simplicity. The MATLAB circuit model for the buck converter can be seen in Figure 3.6. 

The buck converter is a switched mode power supply, meaning it provides digital control 
through PWM of a switch resulting in a linear response in output voltage. It uses a two 
switches, a transistor and a diode, and an inductor to continually maintain a voltage on the 
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Figure 3.5: Bucking DC —DC converter showing simple circuit diagram (a), output voltage 
and frequency spectrum of output (b). Image from [9]. 


output capacitor. Buck converters are commonly used when the input voltage is too high 
and needs to be lowered to an appropriate level. A step down converter produces a lower 
average output voltage than the DC input voltage Vj. Its main application is in regulated 
DC power supplies and DC motor speed control [9]. 

The v 0 i waveform in part (b) of Figure 3.5 is shown as a function of switch position. This 
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Figure 3.6: Overview of the buck converter model in Simulink. 


waveform is what appears at the input to the low-pass filter in part (a) of Figure 3.5. The 
average output voltage can be calculated in terms of switch’s D using 

1 rT s l / rt on rT s \ + 

v 0 = - / Vo(t)dt = -{ v d dt+ 0 dt)= ^ V d = DV d (3.3) 

is Jo is \J 0 Jt on J is 

where, V 0 is the DC output voltage, T s is the switching period, v () (t) is the instantaneous 
output voltage as a function of time, t on is the duration the switch is on per switching period, 
V d is the DC input voltage, and D is the duty cycle of the converter. 

When the active time of the switch t on increases, it increases the DC voltage at the output. 
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The D is related to the ‘on’ time of the switch divided by the period of switching: 


D = 


T s 


(3.4) 


Any linear change in the duty cycle parameter causes a linear change in the output voltage 
[9], 

As stated by Mohan: 


In actual applications, the forgoing circuit has two drawbacks: (1) In practice 
the load would be inductive. Even with a resistive load, there would always 
be certain associated stray inductance. This means that the switch would have 
to absorb the inductive energy and therefore it may be destroyed. (2) Without 
the LC filter, the output voltage fluctuates between zero and Vj, which is not 
acceptable in most applications. [9] 


Normally, the combination of an active inductor serially connected to a capacitor would 
cause arcing across a switch trying to disconnect the source current. To provide an alternate 
current path while the switch has disconnected the source, a diode is included in the circuit. 
By employing the use of a low-pass filter, an inductor and capacitor combination as shown 
in Figure 3.7, the output voltage and current fluctuations delivered to the load are greatly 
reduced. 


The corner frequency f c of this low-pass filter is calculated from 

2kVlC 


(3.5) 
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Figure 3.7: Buck converter in continuous conduciton mode with switch on in part (a) and 
off in part (b). Image from [9] 


The filter cut-off frequency is selected to be much lower than the switching frequency of 
the converter, which reduces the output voltage ripple to a reasonable level [9]. 

During the interval when the switch is ‘on’, the diode is reverse biased, and the input 
provides energy to the load as well as the inductor. During the interval when the switch is 
‘off’, the inductor current flows through the diode transferring some of its stored energy to 
the load. 

This switching can cause some ripple in the output voltage. Due to the nature of inductors, 
the current needs to ramp up until the switching event occurs again. During this time, the 
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small changes in current induces a ripple on the output voltage as seen in Figure 3.8. 
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3.2.3 Load 

Simulating an active load can be very complex. In order to optimize the buck converter and 
MPPT algorithm during simulation, a resistive load that remains constant through each test 
was used. The first order model of the output load can be seen in Figure 3.9. 

The MPPT controller adjusts the D of the buck converter to match the MPP of the array, 
thereby, optimizing output power. By analyzing the output voltage, we can evaluate the 
effectiveness of the buck converter. 
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Figure 3.9: Output load resistor with measurement blocks connected. 


3.2.4 MOSFET 

The MATLAB model as seen in Figure 3.10 for a MOSFET is quite extensive. The ad¬ 
justable parameters are the FET resistance R on , the internal diode inductance L on , the in¬ 
ternal diode resistance R c i, the internal diode forward voltage Vj. the snubber resistance R s , 
and the snubber capacitance C s . 
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Figure 3.10: MOSFET block from MATLAB SimPowerSystems 
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3.2.5 Pulse Width Modulation Generation 


For simulations in MATLAB, a gate drive circuit model is not necessary since the gating 
function assumes an ideal switch. The model for a MOSFET was found in the MATLAB 
blockset SimPowerSystems as seen in Figure 3.10. 

The MOSFET gate control for the model can be seen in Figure 3.11. The percentage of 
time the switch is ‘on’ for a set switching speed is D. We can use the comparison of D with 
a triangle wave that ranges from zero to one. When the value of the triangle waveform is 
higher than D, the MOSFET is ‘off’. Similarly, when the value of the triangle waveform is 
lower than D, the MOSFET is ‘on’. 


Trianglel 



For the MATLAB model, only a binary signal is needed to control the gate. Using this 
simplified idea of the gate, we can compare the duty cycle with a triangle waveform set to 
repeat every 40 /is with the peak of the triangle at 20 fls. This gives a control frequency f s 
of 25 kHz when calculated with 

f*=L < 3 - 6 ) 

When the value of D is less than than the triangle waveform the output goes to ‘one’. 
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3.2.6 MPPT Algorithms 

There are many ways to track an MPP. In this section, two of the most popular MPPT 
algorithms P&O and InC are compared. 

The perturb and observe (P&O) algorithm for finding the MPP is widely used in commer¬ 
cial applications. This method uses current and voltage sensors to measure the output of 
a solar array. The P&O algorithm can be the cheapest to implement. Generally, P&O is 
also the simplest to implement due to its reliance on passive element circuits to measure 
the current and voltage. These measurements are the only outside data sources the MPPT 
algorithm needs to determine whether to increase or decrease the duty cycle. By compar¬ 
ing the input power to the previously measured input power, we can determine whether to 
lower or raise the duty cycle to continue tracking the MPP. 

By changing the duty cycle after a set amount of time, the P&O algorithm checks the 
previously measured input power. If the newly measured input power is greater than the 
previously measured input power and the new input voltage is greater than the old input 
voltage, it decreases D to move closer to the peak power. If the input voltage was less than 
previously measured and it had greater input power, the algorithm increases D. 

If the new power measured is less than previously measured and the new voltage is greater 
than the old voltage, D decreases to converge to an MPP. Finally, if the power and voltage 
are less than previously measured, the algorithm decreases D. This algorithmic flow chart 
can be seen in Figure 3.12. 

There are a few optimizable parameters within this algorithm. First, the f s can affect the 
size requirement for the inductor and capacitor. If f s increases, then the sizes of the inductor 
and capacitor decreases. Secondly, the amount of change in duty cycle can be modulated 
to quickly converge from large tracking errors. This A D is essentially the gain of the con- 
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Figure 3.12: Flow chart of Perturb and Observe Algorithm. Image from [20]. 


troller. If the A D is too high it can overshoot and oscillate indefinitely around the optimal 
point, yet never reach it. 
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Figure 3.13: Flow chart of the InC Algorithm. Image from [21]. 


The incremental conductance algorithm (InC) is another method of MPP tracking. This al¬ 
gorithm uses incremental measurements on the change in conductance of the photovoltaic 
array. By comparing these calculated incremental measurements, we can calculate whether 
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there was an increase or decrease in power since it was last measured. The incremental 
conductance is defined as {Alpy / AVpy). By comparing this measurement to the actual 
conductance of the photovoltaic array, one can determine which side of the MPP the cur¬ 
rent operating point is located. “TnC can track rapidly increasing and decreasing irradiance 
conditions with higher accuracy than P&O. One disadvantage of this is the increased com¬ 
plexity when compared with P&O [20]." 

The algorithm used can be seen as a flow chart in Figure 3.13. This algorithm has fewer 
branches for each node and is less complicated than the P&O algorithm’s flowchart; how¬ 
ever, the simplicity of the algorithm flowchart does not account for the increased complex¬ 
ity of binary calculations required for this algorithm. 

Since calculation of both the InC and conductance requires the use of division, a large 
amount of processing power is necessary in order to acquire the result before the next 
sample time. If the calculation to find the InC carries over into the next sample time, the 
controller will be controlling the buck converter further in the past the longer the system 
runs. If we have each calculation occur as fast as possible by using more power hungry 
processors with the hardware to perform division, we will be able to appropriately use the 
calculations for this algorithm and, thereby, optimize the output power. 
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CHAPTER 4: 

COMPONENT DESIGN 


4.1 Requirements and Description 

For this simulation, the requirements for the proper operation of the buck converter and 
solar array at the MPP are given by a maximum converter input power of 135 W, a converter 
input voltage range of 15 V — 25 V, an output power of 128 IF, and a conversion efficiency 
> 95%. 


4.2 Simulation Data 

By simulating all the designs of the circuit in MATLAB, proof of operation for the buck 
converter acting to track the MPP are produced. Without the ability to simulate the circuit, 
the amount of time required to tune the circuit as well as select the proper components to 
handle the power required would be increased. The procurement procedure for acquiring 
the power components would also incur a great time loss. 

4.3 Calculation of Component Values 

In order to guarantee continuous current mode ( CCM ) for the buck converter, there are a 
few requirements that need to be met in order for the circuit to function as required. 

4.3.1 Input Capacitor Sizing 

In order to store charge while the switch is in the ‘off’ position, a capacitor is placed at the 
input to the buck converter. The capacitor stores charge when the buck switch is in the ‘off’ 
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state, alleviating source current discontinuity. The value calculated using 


r 

'-'it 


(MIN) 


D(l-D) 

Vpp(MAX) x fs 


(4.1) 


where C in ^ MIN ^ is the minimum capacitor size (F), I a is the solar array output current (A), 
Vpp(MAX) is th e maximum voltage from the solar array (V), and f s is the switching fre¬ 
quency (Hz), is 3 flF with a D of 0.5. 


4.3.2 Switching Frequency 

The switching frequency for the MOSFET in the buck converter was computed by inserting 
T s used in simulation, 4 x 10 -5 s, into Equation 3.6. we obtain a switching frequency of 
25 kHz. In order to minimize output ripple, it is essential to have the corner frequency of 
the low-pass filter at the output f c « f s . 


4.3.3 Low-Pass Filter Values 


As a subset of the buck converter, there is an LC low-pass filter present on the output. By 
first solving for L using 


I LB, 


max — 


T s Vd 
8 L ’ 


(4.2) 


with V c i set to 15 V, T s set to 40 /is, and lLB,mca set to 4.1 A, half of the panels’ I sc , we 
get 27 /dH for the inductor. While this value is fine in simulation, there might be a slight 
difficulty finding this specific value in a real-world component. 


From this point we need to size the capacitor to have the filter corner frequency f c be much 
less than f s . Two orders of magnitude in frequency is usually sufficient to remove f s from 
the output. Since f s is 25 kHz, we would like our filter to have a f c of 2000 Hz.. This will 
allow almost pure DC to pass to the output, reducing the ripple in current and voltage. By 
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inputting our previously found value for L into 


fc = 


1 

2 kVlc' 


(4.3) 


we can solve for C. The calculated value for C in the output filter portion of the buck 
converter is 234 \iF. 


4.3.4 Duty Cycle 


The duty cycle D is modulated in our simulation by the controlling algorithm. By varying 
D to control the solar arrays’ operating point, we see that the output voltage is the input 
voltage times D using 


Vo _ ton 

Vd~Y s ' 


(4.4) 


Since the main purpose is to extract as much power from the solar array as possible, it 
is important to note that the output voltage will be unregulated unless energy storage is 
placed at the output of the buck or the resistance is sized such that the buck does not 
become saturated (D = 1). The additional work required will be discussed in Chapter 6. 


4.3.5 Desired Maximum Output Voltage Ripple 

To calculate the maximum voltage ripple of the circuit, we can input the parameters found 
in the previous sections of this chapter into 

AVo= wl (l ~ D)T - <4 ' 5) 

An output ripple voltage of < 1% is desired in most cases. 
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CHAPTER 5: 
RESULTS 


5.1 Validation of MPPT operation 

5.1.1 Irradiance Inputs 

In order to evaluate the performance of each algorithm with a control, three profiles of 
irradiance were used to simulate irradiance on the solar array. The first run was a constant 
irradiance at 1000W /nr. The next run used a step from 600 IT /nr to 1000 W jnr as seen 
in Figure 5.1. The final input to test the systems algorithmic response was a variable level 
of irradiance as seen in Figure 5.2. This final test was to see how well each algorithm 
handled rapidly changing irradiance. 


5.2 Tuning of Duty Cycle Increment. 

5.2.1 P&O 

It is quite easy to tune AD using P&O. Since there is a much longer time between pertur¬ 
bations, approximately 1000 cycles of the simulation, the output power has more time to 
settle before the adjustments are made to the set point of D. This is also aided by the two 
stage gain structure which depends on the error in power. When the difference in power is 
large, then the AD is also large to drive the operating point closer to the MPP. Similarly, 
when there is a small difference in power, then the AD is smaller and overshoots less and 
responds more quickly to small changes in irradiance or temperature. 
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Figure 5.1: The step function used as input to the solar array. 



Figure 5.2: Varying levels of irradiance used as input to the solar array in order to mimic 
passing clouds. 

5.2.2 InC 

While configuring the entire system, the amount of change in D is much smaller than that 
used in the P&O algorithm. Since InC uses samples in time much closer together than 
the P&O algorithm, the AD needs to be smaller to avoid overshooting the target while 
measuring previous changes to the set point. Almost all changes in D will cause a normal 
step response in the output power. By changing the set point of D so frequently, this 
algorithm does have the aspect of being able to compensate for quick changes in shading 
or irradiance. 

Even after reducing the AD to 0.0001, the algorithm still oscillated around the set point 
desired. This oscillation can be attributed to the control loop speed being too fast as well 
as the resolution of AD being too large. The mathematics required in the algorithm has 
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Table 5.1: Average overall converter efficiency for the whole simulation time organized by 
run type. 


Run 

Overall Efficiency 

1 - P&O - Constant 

93.42% 

2 - P&O - Variable 

93.21% 

3 - P&O - Step 

93.42% 

4 - Constant D - Constant 

93.29% 

5 - Constant D - Step 

92.80% 

6 - Constant D - Variable 

92.64% 

7 - InC - Step 

93.68% 

8 - InC - Variable 

93.69% 

9 - InC - Constant 

93.86% 


two checks to see if the change is equal to zero. This can lead to some problems since the 
change between each reading of voltage and current from the output of the solar array can 
sometimes contain noise from either a change in set point of D or from either temperature 
or irradiance change. Further work is required to account for this phenomenon. 


5.3 Simulation Results 


5.3.1 Power Efficiency 


In order to calculate the power efficiency (rj) we need to know the maximum power under 
the given circumstances that the solar array will produce. By integrating both the output 
power measured and the maximum power, we can calculate the MPPT efficiency ( t\mppt ) 
using 


Pmppt 


fo s Pout (t)dt 

j() S PMax(t)dt 


(5.1) 


The results can be seen in Table 5.1. 


Since this is a simulation, there is a great amount of determinism created by running the 
same data with different controls. Essentially, the true MPPT can be extracted from the 
solar array after running the algorithm in order to compute the error in the MPP. 
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Table 5.2: Percentage of output voltage as ripple organized by run type. 


Run 

Output Voltage Ripple % 

1 - P&O - Constant 

0.30% 

2 - P&O - Variable 

0.42% 

3 - P&O - Step 

0.15% 

4 - Constant D - Constant 

0.30% 

5 - Constant D - Step 

0.30% 

6 - Constant D - Variable 

0.25% 

7 - InC - Step 

0.31% 

8 - InC - Variable 

0.23% 

9 - InC - Constant 

0.25% 


5.3.2 Output Voltage Ripple 

The output voltage ripple was calculated from the output voltage versus time plots from 
each run. The key point to realize is that all the algorithms were being controlled by a 
tuned buck converter, and therefore, the output voltage ripple while in steady state was 
well within the specification of AV o /V 0 < 1% as seen in Table 5.2. 


5.4 Comparison of Control Algorithms 

By comparing the output waveforms from each run,we obtain a greater understanding of 
how these algorithms respond to different solar array irradiance changes. 

5.4.1 Run 1: Perturb and Observe with Constant Irradiance. 

This is the first run of the MATLAB model. This run used an input of constant irradiance 
to the solar array. The response curve which shows the output voltage has a measurable 
amount of ripple as seen in Figure 5.4. The current response has a small initial overshoot; 
then the MPPT controlled the PWM of the gate in order to bring it closer to the MPP , 
which can be seen in Figure 5.5. The output power remains almost constant from 0.003 s 
onward as illustrated in Figure 5.6. The converter reaches a steady efficiency of 93.42% 
and also remains almost completely flat for the rest of the simulation seen in Figure 5.3. 
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Figure 5.3: P&O Run 1. Converter efficiency (top) and irradiance (bottom) versus time. 



Figure 5.4: P&O Run 1. Source voltage (top) and output voltage (bottom) versus time. 
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Figure 5.5: P&O Run 1. Source current (top) and output inductor current (bottom) versus 
time. 



Time (s) 



Figure 5.6: P&O Run 1. Source power (top) and output power (bottom) versus time. 
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5.4.2 Run 2: Perturb and Observe with Variable Irradiance. 

This is the second run of the MATLAB model. This run used an input of variable irradiance 
to the solar array. The source and output voltage both track the input irradiance changes 
as seen in Figure 5.8. The current response has a small initial overshoot at each change in 
irradiance. The source current tended to oscillate quite a bit when this controller was in 
operation, which can be seen in Figure 5.5. The output power remains quite close in mag¬ 
nitude compared to the input power, as illustrated in Figure 5.6. The converter efficiency 
shot above 100% because the input power and output power were averaged over 0.0001 s; 
therefore, there was a small lag in the input power. Overall, the efficiency did not drop 
below 80% even with the very abrupt changes in irradiance, seen in Figure 5.3. 
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Figure 5.7: Perturb and Observe Algorithm Run 2. Converter efficiency (top) and irradiance 
(bottom) versus time. 
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Figure 5.8: Perturb and Observe Algorithm Run 2. Source voltage (top) and output voltage 
(bottom) versus time. 
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Figure 5.9: Perturb and Observe Algorithm Run 2. Source current (top) and output inductor 
current (bottom) versus time. 
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Figure 5.10: Perturb and Observe Algorithm Run 2. Source power (top) and output power 
(bottom) versus time. 
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5.4.3 Run 3: Perturb and Observe with Step in Irradiance. 

This is the third run of the MATLAB model. This run has a step in irradiance as input 
to the solar array. The algorithm seems to have corrected in the wrong direction as seen 
in Figure 5.12. The increase in irradiance was supposed to correspond to an increase in 
power output. This transient is most likely due to input noise. The source current tended 
to oscillate only while the irradiance was low, which can be seen in Figure 5.13. The 
output power eventually did recover after around 0.002 s which is the sample rate of the 
algorithm, as illustrated in Figure 5.14. The output power during this incorrect swing was 
around 10 IT. The converter efficiency had a very unstable section that corresponded to 
tracking the incorrect direction of the MPP. The efficiency momentarily dropped to 25%, 
which is illustrated in Figure 5.11. 




Figure 5.11: Perturb and Observe Algorithm Run 3. Converter efficiency (top) and irradi¬ 
ance (bottom) versus time. 
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Figure 5.12: Perturb and Observe Algorithm Run 3. Source voltage (top) and output volt 
age (bottom) versus time. 
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Figure 5.13: Perturb and Observe Algorithm Run 3. Source current (top) and output indue 
tor (bottom) current versus time. 
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Figure 5.14: Perturb and Observe Algorithm Run 3. Source power (top) and output power 
(bottom) versus time. 
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5.4.4 Run 4: Constant D with Constant Irradiance. 

This run was essentially a baseline test for the converter. The source and output voltage 
reach their peak after 0.001 s with very little ripple in the output voltage as seen in Fig¬ 
ure 5.16. The inductor current has a very small amount of ripple, which has a similar 
magnitude as the source current, as can be seen in Figure 5.17. The output power has a 
small dip at 0.001 s due to the resonant nature of the output filter section of the buck con¬ 
verter. The response waveform has a normal damped response as seen in Figure 5.18. The 
efficiency as expected has the same profile as the output power due to the static nature of 
the input irradiance which is illustrated in Figure 5.15. 



Figure 5.15: Constant Run 4. Converter efficiency (top) and irradiance (bottom) versus 
time. 
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Figure 5.16: Constant Run 4. Source voltage (top) and output voltage (bottom) versus time. 
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Figure 5.17: Constant Run 4. Source current (top) and output inductor current (bottom) 
versus time. 


50 






































































Time (s) 



Figure 5.18: Constant Run 4. Source power (top) and output power (bottom) versus time. 
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5.4.5 Run 5: Constant D with Step in Irradiance. 

This run used the constant D controller with a step in the input irradiance. The input 
and output voltage responded predictably to the change in irradiance; however, it was not 
tracking the MPP. The output voltage can be seen in Figure 5.20. The buck converter 
was operating at the MPP while under full irradiance. When the input irradiance reached 
1000 W/m 2 , the output power was 128 W. Since the solar array voltage was so low during 
the first O.OO 85 , the output current also stayed at a lower than normal level, which can be 
seen in Figure 5.21. Since the normal output power at 60% irradiance is 81 W, the panels 
are operating far from the MPP by producing only 50 W, as seen in Figure 5.22. The 
converter efficiency was disrupted by the change in irradiance at 0.008 s, yet the converter 
efficiency remained fairly constant as seen in Figure 5.19. 
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Figure 5.19: Constant Run 5. Converter efficiency (top) and irradiance (bottom) versus 
time. 
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Figure 5.20: Constant Run 5. Source voltage (top) and output voltage (bottom) versus time. 
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Figure 5.21: Constant Run 5. Source current (top) and output inductor current (bottom) 
versus time. 
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Figure 5.22: Constant Run 5. Source power (top) and output power (bottom) versus time. 
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5.4.6 Run 6: Constant D with Variable Irradiance. 

With a variable irradiance, this controller performed the worst of all the algorithms. The 
source and output voltage profiles can be seen in Figure 5.24. At the lowest point in the 
irradiance, the output current dropped to 3.5A. At this time the output voltage was 6 V. 
The output power was approximately 30IT as seen in Figure 5.26. Luckily, the inductor 
current stayed very stable, seen in Figure 5.25, since the D was not changing. The converter 
efficiency was reasonably stable but overall very low as seen in Figure 5.23. 



Time (s) 



Figure 5.23: Constant Run 6. Converter efficiency (top) and irradiance (bottom) versus 
time. 
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Figure 5.24: Constant Run 6. Source voltage (top) and output voltage (bottom) versus time. 
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Figure 5.25: Constant Run 6. Source current (top) and output inductor current (bottom) 
versus time. 
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Figure 5.26: Constant Run 6. Source power (top) and output power (bottom) versus time. 
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5.4.7 Run 7: InC with Step in Irradiance. 

This run used the InC algorithm, and the input irradiance has a step from 600 W / m 2 to 
1000 W/m 2 . While the irradiance is low, the source voltage tracks the MPP quite well and 
keeps the source and output power balanced as seen in Figure 5.28. There was a small 
amount of ripple in the source current as seen in Figure 5.29. Because of the current ripple, 
there is also ripple in the source power measurements. Luckily, the source ripple current 
does not transfer to the output current. The buck converter stabilizes the step in irradiance 
after 0.002 s as seen in Figure 5.30. The converter efficiency does jump around quite a bit 
and has a rather large amount of ripple when it is stabilized as seen in Figure 5.27. This 
oscillation is due to the fast sampling rate set in the algorithm code. 




Figure 5.27: InC Run 7. Converter efficiency (top) and irradiance (bottom) versus time. 
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Figure 5.28: InC Run 7. Source voltage (top) and output voltage (bottom) versus time. 
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Figure 5.29: InC Run 7. Source current (top) and output inductor current (bottom) versus 
time. 
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Figure 5.30: InC Run 7. Source power (top) and output power (bottom) versus time. 
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5.4.8 Run 8: InC with Variable Irradiance. 

This run uses the InC algorithm with a variable input irradiance to the solar array. The 
source voltage and output voltage, seen in Figure 5.32, have a damped response after each 
change in the irradiance. The step size in D can be seen in Figure 5.33, where there is a 
stair step change in the current after 0.008 s. The power output of the converter tracked the 
MPP fairly well as seen in Figure 5.34. The largest oscillation occured in the calculation 
of efficiency. Since it was averaged across 0.001 s there was a false surge in efficiency. 
The inductor and capacitor need time to return to the lower level of output power. This 
also caused the delay in the converter efficiency rising after the irradiance began increasing 
again as seen in Figure 5.31. 
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Figure 5.31: InC Run 8. Converter efficiency (top) and irradiance (bottom) versus time. 


61 































30-1-1-1-1-1-1-1-1-r 



g_i_i_i_i_i_I_i_i_l_ 

0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02 


Time (s) 



Figure 5.32: InC Run 8. Source voltage (top) and output voltage (bottom) versus time. 
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Figure 5.33: InC Run 8. Source current (top) and output inductor current (bottom) versus 
time. 
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Figure 5.34: InC Run 8. Source power (top) and output power (bottom) versus time. 
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5.4.9 Run 9: InC with Constant Irradiance. 

This run uses the InC algorithm with a constant input irradiance . The source voltage took 
longer to reach an acceptable level. The algorithm startes to track the MPP but due to the 
short length of the run, it does not actually reach the MPP. Although, it can be seen in 
Figure 5.36 that the algorithm was tracking in the correct direction and would reach the 
MPP in the next few iterations of the algorithm. The inductor current seen in Figure 5.37 
has four jumps in D while tracking the MPP. With a final output power near 120 W, this 
algorithm takes a few iterations before converging to the MPP as seen in Figure 5.38. The 
converter efficiency is 97% for almost the entire run duration as seen in Figure 5.35. 



Figure 5.35: InC Run 9. Converter efficiency (top) and irradiance (bottom) versus time. 
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Figure 5.36: InC Run 9. Source voltage (top) and output voltage (bottom) versus time. 
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Figure 5.37: InC Run 9. Source current (top) and output inductor current (bottom) versus 
time. 
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Figure 5.38: InC Run 9. Source power (top) and output power (bottom) versus time. 
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CHAPTER 6: 

CONCLUSION AND FUTURE WORK 


6.1 Conclusion 

Overall the comparison of the MPPT algorithms was a worthwhile endeavor. Many of the 
problems inherit in this type of system can be solved or optimized to an acceptable level. 
This work could benefit future development solar powered UAV s. 

In this thesis, the overall operation of the buck converter was covered. The different con¬ 
tributions of each component were discussed and how they each affected the output and 
functionality of the circuit. By optimizing the buck converter independent of the MPPT 
and the solar array, the efficiency of the converter was greater than 96% during steady-state 
operation. This makes the conversion from the solar array quite efficient and attractive. 

The MPPT was discussed after the buck converter and greatly improved the tracking of the 
MPP. By having the ability to track the MPP , the MPPT algorithms further optimized the 
low efficiency solar array output. Since the solar array is the source of power to the system, 
any improvement in the ability to track the MPP , also greatly improves the overall system 
performance. 

Next, the solar array model was discussed as well as designed around a single diode model. 
By modeling the solar array based on measured parameters, the output current and voltage 
were measured and used to drive the buck converter. 

The actual data gathering runs were done in three parts. The first part included the P&O 
algorithm with constant 1000 W/m 2 irradiance, a step in irradiance from 600 W jnr to 
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1000 W/m 1 at 0.08 .s' to 0.01 5, and a random changing of irradiance to the solar array. 
Next, the constant D controller was tested with the same three input signals. The constant 
D was tested as the reference for comparison with the other algorithms. Finally, the InC 
algorithm was tested; although, the tuning of the parameters for this algorithm was the most 
difficult. The algorithm leads to either a very fast response to transients or undershoots the 
MPP. Because of the difficulty optimizing and the amount of time required to actually 
run each test, the InC algorithm was not optimized to match the performance of the P&O 
algorithm. 

By including a two-stage change in I). the P&O algorithm was further optimized. This 
innovation was necessary in order to converge quickly to the MPP yet still have fine gran¬ 
ularity when approaching the MPP. This addition was not very difficult to implement and 
overall provided very fast tracking and convergence to the MPP. In the Appendix, the code 
used for this algorithm is provided. It can be seen that with a multistage MPPT , the overall 
speed of the system will remain fast and track the MPP well. 

Overall, the outcome of the P&O algorithm tests were quite good. Since the algorithm 
can be optimized to meet certain timing requirements, it responded faster and performed 
better than the InC algorithm. This can only be accounted for due to the complexities of 
the InC algorithm computations and speed of operation. Due to the difficulty in optimizing 
the InC algorithm, the results from those runs do not accurately give the optimal operation 
of the algorithm. Given more time to properly optimize the AD, the algorithm should have 
tracked the MPP quicker and with more certainty than the P&O algorithm. 

6.2 Future Work 

There are a few optimizations that were not integrated into the thesis but would be valuable 
to include. 
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6.2.1 UAV Flight Endurance 

In the modern UAV, most onboard components are run from a battery. Similar to satellites, 
UAV s require battery power to be replenished at some point. If all of the power for the sys¬ 
tem comes from the battery, it seems logical to try and extend the range and time the UAV 
will be on a mission. In order to accomplish this, thin film solar panels can be laminated to 
the surface of the wings of these UAV s and used to deliver power to the batteries while in 
flight. With the inclusion of an MPPT , the solar panels on the wing will be optimized for 
when the UAV turns to a non-optimal power generation angle. 


6.2.2 Second Stage Constant Voltage Output 

In order to actually use the output power from the MPPT , a second stage constant voltage 
controller would also be necessary to provide the power required to charge a battery system. 
This battery system would then have all the equipment connected to that so the voltage 
would be quite stable at the equipment. 

Without this extra controller the battery charger attached to the buck converter would not 
have a high enough voltage if the irradiance dropped too low. Depending on the amount 
of light and the configuration of the solar panels, a buck-boost converter might be better 
suited for this application. 


6.2.3 Test MPPT on Real Hardware 

A proof-of-concept for this simulation would be to use real hardware to control a buck 
converter. Unfortunately, we did not have the required time to change the laboratory setup 
and run the algorithm on real hardware. 
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6.2.4 More Complex Variable AD Modulation 

Even more advanced than the two stage P&O D modulation would be to have the A D be 
variable based on a lookup table of the MPP. This would have an adjustable gain which 
was adjusted by the distance to the MPP. The two-stage P&O algorithm implemented in 
this thesis is quite similar to automatic sliding AD modulation except that it does not adapt 
to changing parameters as well. 


70 



Appendix: DATASHEET & MATLAB CODE 

Datasheet for Kyocera KD135GX-LP 
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The New Value Frontier 


v 


MODEL 

KD135GX-LP 



KD135GX-LP 


HIGH EFFICIENCY MULTICRYSTAL 
PHOTOVOLTAIC MODULE 

© 

LISTED 

HIGHLIGHTS OF 

KYOCERA PHOTOVOLTAIC MODULES 

Kyocera’s advanced cell processing technology 

and automated production facilities produce a highly efficient 

multicrystal photovoltaic module. 

The conversion efficiency of the Kyocera solar cell is over 16%. 
These cells are encapsulated between a tempered glass cover 
and a pottant with back sheet to provide efficient protection 
from the severest environmental conditions. 

The entire laminate is installed in an anodized aluminum frame 
Equipped with plug-in connectors. 


3 KyocERa 



to provide structural strength and ease of installation. 


APPLICATIONS _ 

KD135GX-LP is ideal for grid tie system applications. 

• Residential roof top systems • Water Pumping systems 

• Large commercial grid tie systems • High Voltage stand alone systems 

• etc. 


QUALIFICATIONS 

• MODULE : UL1703 listed • FACTORY : ISO9001 and ISO 14001 

QUALITY ASSURANCE _ 

Kyocera multicrystal photovoltaic modules have passed the following tests. 

• Thermal cycling test • Thermal shock test • Thermal / Freezing and high humidity cycling test • Electrical isolation test 

• Hail impact test • Mechanical, wind and twist loading test • Salt mist test • Light and water-exposure test • Field exposure test 


LIMITED WARRANTY 

1 year limited warranty on material and workmanship 
r&20 years limited warranty on power output: For detail, please refer to "category IV" in Warranty issued by Kyocera 

(Long term output warranty shall warrant if PV Module(s) exhibits power output of less than 90% of the original minimum rated power specified at the time of sale within 
10 years and less than 80% within 20 years after the date of sale to the Customer. The power output values shall be those measured under Kyocera’s standard 
measurement conditions. Regarding the warranty conditions in detail, please refer to Warranty issued by Kyocera) 


ELECTRICAL CHARACTERISTICS 


Current-Voltage characteristics of Photovoltaic 
Module KD135GX-LP at various cell temperatures 



Current-Voltage characteristics of Photovoltaic 
Module KD135GX-LP at various irradiance levels 
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SPECIFICATIONS 


KD135GX-LP 


■ Physical Specifications 


Unit: mm (in.) 


668 (26.3in.) 



35.7 
(1.41 in.) 



36(1.42in.) 


■ Specifications 


■ Electrical Performance under Standard Test Conditions (*STC) 

Maximum Power (Pmax) 

135W (+5%/-5%) 

Maximum Power Voltage (Vmpp) 

17.7V 

Maximum Power Current (impp) 

7.63A 

Open Circuit Voltage (Voc) 

22.1V 

Short Circuit Current (isc) 

8.37A 

Max System Voltage 

600V 

Temperature Coefficient of Voc 

-.080 V/”C 

Temperature Coefficient of Isc 

5.02x10~ 3 A/°C 


*STC : Irradiance 1000W/m 2 , AMI.5 spectrum, cell temperature 25‘C 


■ Electrical Performance at 800W/m 2 , *N0CT, AMI.5 

Maximum Power (Pmax) 

95W 

Maximum Power Voltage (Vmpp) 

15.6V 

Maximum Power Current (impp) 

6.10A 

Open Circuit Voltage (Voc) 

19.9V 

Short Circuit Current (isc) 

6.82A 


*NOCT (Nominal Operating Cell Temperature): 49'C 


■ Cells _ 

Number per Module 36 


■ Module Characteristics 


Length X Width X Depth 

1500mm(59.1in)x668mm(26.3in)x36mm(1,4in) 

Weight 

13.0kg(28.7lbs.) 

Cable 

(+)760mm(29.9in),(-)1840mm(72.4in) 



■ Junction Box Characteristics 


Length X Width x Depth 

100mm(3.9in)x 108mm(4.3in)x 15mm(0.6in) 

IP Code 

IP65 



■ Others 


"Operating Temperature 

-40"C~90"C 

Maximum Fuse 

15A 

‘This temperature is based on cell temperature. 



Please contact our office for further information 


^KyocERa 

KYOCERA Corporation 


■ KYOCERA Corporation Headquarters 

CORPORATE SOLAR ENERGY DIVISION 
6 Takeda Tobadono-cho 
Fushimi-ku, Kyoto 
612-8501, Japan 

TEL:(81 >75-604-3476 FAX:(81 >75-604-3475 
http://www.kyocera.com/ 

• KYOCERA Solar, Inc. 

7812 East Acoma Drive 
Scottsdale, AZ 85260, USA 

TEL:(1 >480-948-8003 or (800)223-9580 FAX:(1 >480-483-6431 
http://www.kyocerasolar.com/ 

• KYOCERA Solar do Brasil Ltda. 

Av. Guignard 661, Loja A 

22790-200, Recreio dos Bandeirantes, Rio de Janeiro, Brazil 
TEL:(55)21 -2437-8525 FAX:(55)21 -2437-2338 
http://www. kyocerasolar.com.br/ 

• KYOCERA Solar Pty Ltd. 

Level 3, 6-10Talavera Road, North Ryde 
N.S.W. 2113, Australia 

TEL:(61 >2-9870-3948 FAX:(61 >2-9888-9588 
http://www. kyocerasolar.com.au/ 

• KYOCERA Fineceramics GmbH 

Fritz-Muller-Strasse 107, 73730 Esslingen Germany 
TEL:(49)711 -93934-999 FAX:(49)711 -93934-950 
http://www. kyocerasolar.de/ 
solar@ kyocera.de 


• KYOCERA Asia Pacific Pte. Ltd. 

298 Tiong Bahru Road, #13-03/05 
Central Plaza, Singapore 168730 
TEL:(65)6271 -0500 FAX:(65)6271 -0600 

• Kyocera Asia Pacific Ltd. 

Room 801-802, Tower 1, South Seas Centre, 

75 Mody Road, Tsimshatsui East, Kowloon, Hong Kong 
TEL:(852)2723-7183 FAX:(852)2724-4501 

• KYOCERA Asia Pacific Pte. Ltd., Taipei Office 

10F, No. 66, Nanking West Road, Taipei, Taiwan 
TEL:(886) 2-2555-3609 FAX:(886)2-2559-4131 

• KYOCERA (Tianjin) Sales & Trading Corp. 

(Beijing Office)Room 2107, Beijing Huabin International Building, 
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’P&O’ 


function Duty = MPPT(V,I) 

persistent PowerNew PowerOld DifferenceOfPower d DutyDifference n i OldV; 

if isempty(V) 

V=20; 

end 

if isempty(I) 

1 = 0 ; 

end 

if isempty(OldV) 

01dV=0; 

end 

if isempty(PowerOld) 

Power01d=0; 

end 

if isempty(PowerNew) 

PowerNew=0; 

end 

if isempty(DifferenceOfPower) 

DifferenceOfPower=0; 

end 

if isempty(d) 
d=l; 

end 

if isempty(DutyDifference) 
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DutyDifference=0; 


end 

if isempty(n) 
n=500; 

end 

if isempty(i) 

i=l; 

end 

PowerOld = PowerNew; 

PowerNew=V*I; 

DifferenceOfPower = PowerNew-PowerOld; 
if ( i > n ) 
i = 1; 

if (DifferenceOfPower > .01) 
if (V > OldV) 

if(DifferenceOfPower > 1) %if power>l 
DutyDifference=.75; 
d=d+DutyDifference; 
else %if .01<power<l 

DutyDifference=.025; 
d=d+DutyDifference; 

end 

else %V<01dV 

if(DifferenceOfPower > 1) 

DutyDifference=-.75; 
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d=d+DutyDifference; 
else 7„if .01<power<l 

DutyDifference=-.025; 
d=d+DutyDifference; 

end 

end 

else 

if(DifferenceOfPower < -.01) 
if(V<01dV) 

if(DifferenceOfPower <-l) 
DutyDifference=-.75; 
d=d+DutyDifference; 
else °/„if -,01<power <-l 
DutyDifference=-.025; 
d=d+DutyDifference; 

end 

else 7„V>01dV 

if(DifferenceOfPower < -1) 
DutyDifference=.75; 
d=d+DutyDifference; 

else 

DutyDifference=.025; 
d=d+DutyDifference; 

end 

end 

else 
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DutyDifference=0; 


end 

end 


else 

i=i+l; 

end 

Duty=d/(d+1); 

%keep spikes in duty cycle under control. 
if (Duty >.9) 

Duty=.9; 

else 

if(Duty < .1) 

Duty = .1; 

end 

end 

01dV=V; 
end 


’InC’ 

function Duty = MPPTIC(V.I) 

persistent DeltaVoltage DeltaCurrent d DutyDifference n i OldV Oldl; 
if isempty(¥) 
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V=20; 


end 

if isempty(I) 
1 = 0 ; 


end 

if isempty(OldV) 
01dV=0; 


end 

if isempty(Oldl) 
01dl=0; 


end 

if isempty(DeltaVoltage) 
DeltaVoltage=0; 


end 

if isempty(DeltaCurrent) 
DeltaCurrent=0; 


end 


if isempty(d) 
d=l; 


end 

if isempty(DutyDifference) 
DutyDifference=0; 


end 

if isempty(n) 
n=l; 
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end 


if isempty(i) 

i=l; 

end 

DeltaVoltage = V-OldV; 

DeltaCurrent = I-01dI; 

%start main loop 
if ( i > n ) 

i = i; 

if (abs(DeltaVoltage) <= .01) 

if(abs(DeltaCurrent) <= .01) 

°/„No Change; 

else 

if(DeltaCurrent > 0) 
d=d+.0001; 

else 

d=d-.0001; 

end 

end 

else 

if(abs((DeltaCurrent/DeltaVoltage)+I/V) < .001) 
7„no change 

else 

if((DeltaCurrent/DeltaVoltage) < -I/V) 
d=d+.0001; 
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else 


d=d-.0001; 

end 

end 

end 

else 

i=i+l; 

end 

Duty=d/(d+1); 
if (Duty >.9) 

Duty=.9; 

else 

if(Duty < .1) 

Duty = .1; 

end 

end 

01dV=V; 

01dI=I; 

end 
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